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Abstract 

We show that a closed piecewise-linear hypersurface immersed in IR n (n > 3) is 
the boundary of a convex body if and only if every point in the interior of each 
(n — 3)-face has a neighborhood that lies on the boundary of some convex body; no 
assumptions about the hypersurface's topology are needed. We derive this criterion 
from our generalization of Van Heijenoort's (1952) theorem on locally convex hy- 
persurfaces in W 1 to spherical spaces. We also give an easy-to-implement convexity 
testing algorithm, which is based on our criterion. For M 3 the number of arithmetic 
operations used by the algorithm is at most linear in the number of vertices, while 
in general it is at most linear in the number of incidences between the (n — 2)-faces 
and (n — 3)-faces. When the dimension n is not fixed and only ring arithmetic is 
allowed, the algorithm still remains polynomial. Our method works in more general 
situations than the convexity verification algorithms developed by Mehlhorn et al. 
(1996) and Devillers et al. (1998) - for example, our method does not require the in- 
put surface to be orientable, nor it requires the input data to include normal vectors 
to the facets that are oriented "in a coherent way". For R 3 the complexity of our 
algorithm is the same as that of previous algorithms; for higher dimensions there 
seems to be no clear winner, but our approach is the only one that easily handles 
inputs in which the facet normals are not known to be coherently oriented or are 
not given at all. Furthermore, our method can be extended to piecewise-polynomial 
surfaces of small degree. 
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1 Introduction 



Blum and Kannan (1989) suggested a paradigm of output verification. Since a complete 
check of a program is often difficult or not possible - for example, when the source code 
has not been made public - it is important to have algorithms that verify key properties of 
mathematical objects generated by programs. Instead of the source code verification one can 
try to verify the properties of the output that are deemed essential by users of the program. In 
computational geometry this paradigm was developed, among others, by Melhorn et al. (1996, 

1999) and Devillers et al. (1998). For example, the LEDA C++ library contains programs 
verifying convexity of a polygon, Delaunay property of a tiling, etc (Mehlhorn and Naher, 

2000) . Devillers et al. (1998) argue that it is easier to evaluate the quality of the output of 
a geometric algorithm, than the correctness of the algorithm or program producing it. This 
paper contributes to the problem of verification of convexity of a large class of piecewise- 
linear (PL) hypersurfaces in M. n for n > 3. The novelty of our approach is in reducing the 
verification of global convexity of a PL-hypersurface to the verification of local convexity at 
the faces of small codimension. We show that a closed bounded PL-hypersurface realized in 
M. n (n > 3) without local self-intersections is the boundary of a convex body if and only if 
each (n — 3)-face of the hypersurface has a point in its relative interior such that a small 
Euclidean ball at this point is cut by the hypersurface into two pieces, one of which is convex. 
The local convexity condition can also be expressed as that the point has a neighborhood 
on the hypersurface, which lies on the boundary of a convex body. If the hypersurface is not 
bounded, this criterion is invalid. However, if it is known that the hypersurface has at least 
one point of strict convexity - i.e., a point s such that a small ball, centered at s, intersects 
the hypersurface over a set S, which lies on the boundary of a convex body and such that 
S"\s lies in an open half space with respect to some hyperplane through s - our test can still 
be used. In fact, our result is slightly more general: we prove our criterion for any cellwise- 
flat PL-realization in M. n (n > 3) of a semiregular CW-partition of an (n — l)-manifold. The 
technical terms used in this formulation are defined in the next section. 

In this paper we also construct an algorithm for convexity testing that can be applied to any 
closed PL-hypersurface. Our approach does not require any preliminary knowledge of the 
topology of the surface, or any information about its orientability. The direct comparison of 
the complexity of our algorithm and those of Mehlhorn et al. (1996b; 1999) and Devillers 
et al. (1998) is not quite meaningful, since these authors make the following simplifying 
assumptions. 

(CO) The input is known to be an orientable closed hypersurface. The normals to the (n — 1)- 
faces are given as part of the input, and they are all oriented either outwards or inwards 
("coherent orientation") 
(S) the cell-partition is simplicial. 

The last condition is not needed for these algorithms to work; however, it does affect the 
complexity analysis, which is given in the papers only for the simplicial case. It is not 
clear to us why the assumption (CO) is natural; in any case, they significantly simplify all 
considerations. Our approach does not require the assumption (CO). For M 3 , if conditions 
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(CO) and (C) are met, both our and previous approaches have the same complexity, which 
is big-0 in the number of vertices. While Devillers et al.(1998) claimed 0(/o) as the running 
time of for any dimension, this was clearly a typo: first, this bound is impossible, even for 
PL-spheres, by fundamental counting theorems of polyhedral combinatorics and second, the 
pseudocode in (Devillers et al, 1998: Section 3) has running time, in the notation of that 

d-2 d 

paper, of £fj + £ T,fj-i,j(F). 

3=0 {FeT I dimF=d-2}i=3 

Formal definitions and notation are given in Section [2J The paper presu pposes some familiar- 



ity w ith partially ordered sets (pos ets) , linear algebra and geometry fe.g. lKostrikin and Maninl . 



1988 ). basic convexity theory (e.g. Rockafellar . 1990), genera l topol ogy (e.g. Dugundji . 19651 ). 



and basic combinatorial topology (e.g. ISeifert and Threlfalll . Il980l ). 



1.1 On efficiency of Convexity Checkers 

Suppose all computations are done with floating point arithmetic. More formally, consider 
the random access machine (RAM) with unit cost model of computation, where all four 
arithemtic operations are included in the instruction set. As usual, f; j denotes the number of 
incidences between z-dimensional and j-dimensional faces and fj = fj , stands for the number 
of i-dimensional faces. Con sider the case o f R 3 . I f condition (CO) holds, the algorithms by 
Mehlhorn et~aT1 fll996h and iDevillers et all (1l998h . have the same time-complexity of O(fo). 



Our algorithm also has the complexity of O(f ), without requiring the assumption (CO); 
furthermore, its working does not depend on whether the surface is orientable or not. In 
a more general situation, where the conditions (CO) and (S), or some of them, cannot be 
assumed, the complexity of our algorithm, as well as the algorithms of the previous authors 
heavily depend on the following three factors: i) the combinatorics of the cell-partition (e.g. 
simplicial or not), ii) the geometry of the realization (e.g. generic positions of the vertices 
vs. completely general case), iii) the form of the input. 

Regarding iii), for example, the combinatorial information about the input can be given in 
the form of the complete poset of faces, or some subposet of faces, such as, e.g. the vertex- 
facet graph. Furthermore, certain additional topological information (e.g. the knowledge that 
the hypersurface is orientable, or a cyclic order of the facets at each (n — 3)-face) might speed 
up the convexity verification. The geometry of the realization can be given by the equations 
of the facets, or by "coherent" inequalities for the facets (CO), or by positions of the vertices, 
or in the form of inner normals at (n — 3)-cells to the (n — 2)- and (n — 3)-cells. 

Now let us consider the general problem of convexity verification. At one end is the simplified 
setup, where the input hypersurface is simplicial and the realization is sufficiently generic 
so that floating point arithmetic can safely be used. Under these assumptions everything is 
fast, regardless of the method used. On the other end is the completely general setup, where 
nothing can be assumed. One can also consider "intermediate" models, such as, e.g., where 
the input hypersurface is simplicial, but the positions of the vertices are not necessarily 
generic. Another reasonable assumption would be that the hypersurface is not necessarily 
simplicial, but the realizations of the facets are known to be convex. We give our algorithm 
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for the most general case, where nothing can be assumed. One of the motivations for this 
generality is the work of Joswig and Ziegler (2004), who clearly demonstrated that from the 
complexity theory point of view the convex hull problem is most interesting when we cannot 
assume that the vertices (or hyperplanes) are in general position - or that the dimension of 
the space is a small fixed number. 

If sufficient linear-algebraic data and face incidence information are given about the stars of 
(n — 3)-faces (see Section [6731 for details), then the complexity of our algorithmic approach is 
still polynomial (in the Turing machine model). The distinguishing features of our approach 
are: 1) locality of testing and 2) independence from the global topology of the input surface. 
Moreover, our approach generalizes to piecewise-algebraic surfaces (this work is in progress). 



2 Definitions and Background 



From now on X n is used to denote IR n or S n . By a subspace of X n we mean an affine subspace 
in the case of M. n , and the intersection of S n C R n+1 with a linear subspace of M n+1 in the 
case of S n . A hyperplane is a subspace of codimension one. A set K C X n is called convex 
if for any x,y G K there is a minimal geodesic segment with end-points x and y that lies in 
K. Then, the dimension of K, dim K, is understood as the dimension of a minimal subspace 
containing K (since such a subspace is unique for both W 1 and S n , we denote it by aff K 
even in the spherical case). We use int K to denote the interior of K in aff K: in other words 
int K stands for the interior of K relative to aff K. A hypersurface in X n is a pair (^#, r 
where is a manifold of dimension n — 1, with or without boundary, and r : jtft — > X 
is a continuous realization map. Unless specified otherwise, a manifold will always mean a 
manifold without boundary. When the realization map is fixed we do not mention it every 
time. A realization r is called complete if any sequence on which is Cauchy with respect 
to the r-metric on (see Section [3]), converges to a point of ' . A map i : — > X™ 
is called an immersion if i is a local homeomorphism, in such a case we may also refer 
to i) as an immersed hypersurface. A map / : ^ — > X n is called an embedding if / 
is a homeomorphism onto /(*/#). Obviously, an embedding is an immersion, but not vice 
versa. We call a £;-submanifold (with or without boundary) S of ^ flat with respect to the 
realization map r if r : S — > r(S) is an embedding into a A;-subspace of X n . A submanifold is 
called open if it is non-compact and without boundary. 

A convex body in X™ is a closed convex set of full dimension. The hypersurface (^#, r) is 
called locally convex at p e j$ if p has a neighborhood JV P C ~# such that r : jY v — > r(^AQ 
is a homeomorphism and r( t A / p ) lies on the boundary of a convex body K p = K p (jV p ). Often, 
when it is clear from the context that we are discussing the properties of r near p = r(p), 
we say that r is convex at p. If K p can be chosen so that K p \ r(p) lies in an open half-space 
defined by some hyperplane passing through r(p), the realization r is called strictly convex 
at p. We will als o sometimes refer to (^#, r) as strictly convex at r(p). By a theorem of 



Busemannl (119581 ) JV V and K P (^V P ) can always be chosen to allow a support hyperplane H 
at r(p) such that the orthogonal projection of r{^Y p ) onto H is an open (n — l)-ball. Due to 
the local nature of this theorem, it holds in all spaces of constant curvature. When K p and 
jV p satisfy this assumption, we refer to K p as a convex witness for p. 
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Let us recall (e.g. Rockafellar . 1990h that a point p on the boundary of a convex set K is 
called exposed if K has a support hyperplane that intersects K, the closure of K, only at p; 
also, p is called extreme if it does not belong to the interior of any interval contained in dK. 
Thus, an exposed point on a convex body B is a point of strict convexity on the hypersurface 
dB. Conversly, a point of strict convexity p £ ^# for (^#, r) is an exposed point for a convex 



witness K p . Local convexity can be defined in many other, non-equivalent, ways (e.g., see 
van Heijenoort, 1952). 



We will say that the hypersurface (^#, r) is the boundary of a convex body K C X n if r is a 
homeomorphism from ^# onto dK. Hence, we exclude the cases when r(^#) coincides with 
the boundary of a convex body but r is not injective. Our algorithm for PL-hypersurfaces will 
always detect a violation of the immersion property; in the case where r(^) is the boundary 
of a convex body, but r is not a homeomorphism, it will produce the negative answer without 
trying to determine if r(^) is the boundary of a convex body. Of course, the algorithmic 
and topological aspects of this case may be interesting to certain areas of geometry, such as 
origami. Note that for n > 3 a closed (n — l)-manifold j$ cannot be immersed into M. n by a 
non-injective map r so that r(^#) is the boundary of a convex set - any convex hypersurface 
in MJ 1 is simply-connected and any covering map onto a simply-connected manifold must be 
a homeomorphism. However, such immersions cannot be easily ruled out in the hyperbolic 
space HP, as there are infinitely many topological types of convex hypersurfaces in HP for 
n > 2 (Kuzminykh, 2005). 

This paper is mainly focused on convexity of piecewise-linear (PL) hypersurfaces, in particu- 
lar, boundaries of polytopes. Denote by B d the closed unit ball at the origin in M, d . A (disjoint) 
countable partition 7 of a topological space j$ is called a semiregular cell-partition if (1) 
each element C £ 7, called a cell of 7), is homeomorphic to int B dimC , where dim C £ N 
and dimC < dim^#; (2) the closure C (in ^) of each C £ 7 is the union of C and cells 
of smaller dimensions; (3) for each C £ 7 there is a mapping rc '■ C — > B dimC which is a 
homeomorphism onto rc(C) and whose restriction to C is a homeomorphism onto int B dimC . 
Authors that prefer to deal with closed cells refer to cell-partitions as cell- complexes. 



If each cell is contained in the closure of finitely many cells, the partition is called star- 
finite. StarF denotes the subcomplex that consist of all (relatively open) cells whose closure 
contains F. If the closure of each cell is the union of finitely many cells, the partition is 
called closure-finite. When a partition is both closure- and star- finite, it is called locally- 
finite. Often, in the definition of cell-partition (cell-complex) one insists on that each closed 
cell is the image of a closed ball, which forces the compactness for each closed cell - we 
do not make such a requirement. Hence, in our definition the closures of the cells can be 
"semiclosed-semiopen" . Our notion of semiregular cell-partition is a natural generalization of 
the standard notion of regula r cell-partition , also known as regular CW-complex, introduced 
by J.H.C. Whitehead (see e.g. IZieglerl . |2002| ) . Namely, a regular cell-partition is a semiregular 
locally-finite cell-partition, where the closure of each cell is homeomorphic to a closed ball. 



According to our definition, for example, the vertical projection on the plane of the graph (in 
IR 3 ) of a continuous piecewise-linear (more formally, piecewise-affine) function / on IR 2 , which 
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is defined by finitely many affine equations and inequalities, naturally induces a semiregular 
cell-partition CP/ of R 2 : each 2-cell is a maximal linearity set of the function. However, CPy fails 
to be a regular cell-partition due to its unbounded cells. Semiregular partitions are especially 
well-suited for the study of topology of real semialgebraic and subanalytic sets - any such 
set has a canonical finite semiregular cell-partition. 

It is easy to see that any regular cell-partition can be subdivided into a triangulation. Thus, 
a manifold admitting a regular cell-partition belongs to the category of PL-manifolds. It 
can be proven that a semi-regular locally-finite cell-partition can also be subdivided into 
a triangulation. However, since the algorithmic part of this paper deals only with finite 
partitions, we omit this theorem. 

With any semiregular cell-partition there is a natural structure of poset. Namely, for cells F 
and C we have F ■< C if and only if F C C; we write F -< C if F C C but F ^ C. If F ■< C, 
we say that F is a face of C. We will use the same symbol CP for a partition and its poset. It 
is convenient to augment the poset CP with an infinum _L, which is 0, and a supremum T, 
which can be thought of as all of . We use Skd(^#, CP) to denote the (/-skeleton of CP) 
- i.e. the subcomplex of (^#, CP) that consists of all cell of dimensions not exceeding d. 

A subset of R ra is called polyhedral if it is defined by a propositional formula in the language 
of the reals (R) that uses only affine equations and inequalities. A subset S of S n C R ra+1 is 
called polyhedral if S = S n D E, where E C R n+1 is defined by a propositional formula in the 
language of the reals (R) that uses only linear (homogeneous) equations and inequalities. A 
PL-hypersurface in X n is a triple (^#, CP; r), where ^ is a manifold with a semiregular cell- 
partition CP, and r : ^# — > X n is a continuous realization map, such that for each C € CP the 
set r(C) is polyhedral and is homeomorphic to C. Note that although r need not even be an 
immersion, the restriction of r to the closure of any cell C of ^# must be an embedding. The 
realization is called cellwise-flat if the dimension of the affine span of r(C) is equal to dimC. 
Although the term face is used both for abstract cells and their geometric realizations, we 
usually apply it for the realizations. If dimC = k, then r(C) is called a k-f&ce of (~-#,CP;r). 
Throughout the paper all faces, just as all cells in topological partitions, are assumed to 
be relatively open. At times we refer to (n — l)-faces as facets, (n — 2)-faces as ridges, and 
(n — 3)-faces as corners. We may also use these geometric names for the underlying cells of 
CP) - the meaning will always be clear from the context. 

Let us consider a connected PL-hypersurface (^#, CP; r), where r : ^# — > R n maps each cell 
onto a set of the same affine dimension, i.e., dimC = dimaff r(C). Suppose r : ^# — > R n has 
at least one point of strict convexity; also, suppose that r is locally convex at all points of all 
corners of (^, CP). Notice that if the last condition holds for some point of a corner, it holds 
for all points of the corner. Our main Theorem [H] states that under these conditions r is a 
homeomorphism on the boundary of a convex body. This theorem implies a test for global 
convexity of a PL-hypersurface that proceeds by checking the local convexity on each of the 
corners. The pseudo-code for the algorithm is given in Section 16.31 The complexity of this 
test depends not only on the model of computation, but also on the way the surface is given 
as input data. Let the input be the poset of faces of dimensions n — 1 (facets) n — 2 (ridges) 
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and n — 3 (corners). Suppose for each corner-ridge incidence (C, R) we are given a Euclidean 
inner normal to r(R) at r(C), and for each corner-facet incidence (C,F) we are given a 
Euclidean inner normal to r(F) at r(C). If we adopt the algebraic complexity model where 
each of scalar operations {comparison, addition, subtraction, multiplication, and division} 
has unit cost, the complexity of the algorithm is 0(nf n _ 3n _ 2 ) = 0(nf n _ 3n _i). Complexity 
under other models is studied in Section [TJ we also study the complexity of extracting the 
required input information from more common input representations. 

The algorithm consists of f n _3 independent subroutines corresponding to the (n — 3)-faces, 
each with complexity not exceeding big-0 in the number of (n — 2)-cells incident to the 
(n — 3)-face. In addition to the algorithmic implications, our generalization implies that any 
(n — 3)-simple compact PL-hypersurface in M n is the boundary of a convex polytope. 



3 Geometry of Locally-Convex Immersions 

Recall that a path joining points x and y in a topological space T is a continuous map 
p : [0, 1] — > 7. Such a path is called an arc if p is injective. Denote by Arcs ^(x, y) the set of 
all arcs joining x, y G j& . 

An immersion i : ^# — > X n induces a metric di on <M by 

di(x,y) = inf \i(a)\ 

a&Aics^(x,y) 

where \i(a)\ G RUoo stands for the length of the z-image of an arc a joining x and y on ' . 
This metric is called the z'-metric. Of course, for a general continuous realization r it is not 
clear a priori that there is a path of finite length on r(^#) joining r(x) and r(y). That is 
why we need Lemma [H 

Lemma 1 If i is an immersion, then any two points of j$ can be connected by an arc of a 
finite length. In particular, j$ is not only connected, but also arcwise connected. 

The following lemma is used implicitely throughout our proofs and is also important for 
understanding Van Heijenoort and Jonker & Norman' arguments that we are employing. 

Lemma 2 Ifiis an immersion, then the metric topology defined by the i-metric is equivalent 
to the original topology on 

Van Heijenoort's (1952) proofs of these two lemmas, given for X n = R n , work for X n = S n 
without changes. Note that since i is an immersion, then for a "sufficiently small" subset S 
of ^# the map i\§ is a homeomorphism and, therefore, the topology on S that is induced 
by the metric topology of X n , is equivalent to the intrinsic topology of S, i.e., the subspace 
topology. Thus, for sufficiently small subsets of ^ the three topologies considered in this 
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section are equivalent. This fact will be used on numerous occasions without an explicit 
reference to the above lemmas. We will also need the following theorem. 



Theorem 3 (Van Heijenoort, 1952) If a complete locally convex immersion f of a connected 
manifold ^# (dim jtf. = n — 1 ) into M. n (n>3) has a point of strict convexity, then f is a 
homeomorphism onto the boundary of a convex body. 

Recall that a subset C of the linear space M n+1 is called a cone if AC C C for any A G M + . 
A cone C is called pointed if C C. A cone is called salient if it does not contain any linear 
subspace except for 0. If S C M n+1 , then we denote by p ■ S the cone with apex p over S. 

Let x be a point on ^# and let S be a subset of X n such that xDr~ 1 (S) ^ for a realization 
map r : J% — ► X n . We denote by r~ l S the connected component of r~ 1 (S) that contains x. 
Consider now the case of X n = S n . We assume that § n is embedded as the standard unit 
sphere into M n+1 . For a point a G S n C M n+1 we denote by c a the central projection mapping 
from the half of § n that is centered at a onto the tangent plane T a C M n+1 to S n at a. To 
simplify the visual appearence of formulas we will use c x instead of c r ^ in a context where 
the map r : jft. — > § n is fixed. 

The following theorem, whose proof follows the approach taken by Jonker & Norman (1972), 
shows that for spherical immersions absence of a point of a strict convexity cannot result in 
the loss of global convexity, as it happens in the Euclidean case. 

Theorem 4 Let i : ^# — > § n (n >3) be a locally convex complete immersion of a connected 
(n — 1) -manifold ^ . Then (^#, %) is strictly locally convex in at least one point, or = 
d(S n H C), where C is a pointed convex cone in R n+1 . 

In the latter case we say that i is conical. 

Proof. If j$ has a point p such that i(p) is an e xtreme point for so me convex witness K p , 
then, by a classical theorem of Straszewicz (e.g. Rockafellar . 1990h . either i(p) is also an 



exposed point, or there are infinitely many exposed points of K p arbitrarily close to i(p). 
Since an exposed point is a point of strict convexity, the surface i) is strictly convex in 
at least one point and the theorem follows. 

Otherwise, suppose that for all p G ~# the image i(p) is an interior point of a segment on 
Below we will show that this assumption implies that i is a homeomorphism onto 
d(S n H C) for some pointed convex cone C. If a segment through i(jp) cannot be extended 
to a circle on i(^), which is the z-image of a closed curve through p on then pick an 
end point of a maximal segment I p through i(p) and call it i(a). Since i(a) is not extreme, 
it must lie in the interior of another segment. There is a support plane H a for a such 
that int I p and H a do not intersect near i(a). Consider now i~ l H a , which is either a closed 
submanifold of dimension at least 1 or a closed submanifold with boundary. If i^ 1 H a is a 
submanifold with boundary, then the central projection c a i{i~ l H a ) of i{i~ l Ha) on T a = W 1 
must have parallel lines on the boundary, which correspond to intersecting half-circles on § n . 
Furthermore, through any two points on the relative boundary of c i(z~ 1 i7 ) there are two 
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parallel lines contained in the relative boundary of Caify^Ha). Let R C T a be a maximal 
subspace contained in c a i{i~ l H a ). The c a -preimages of points at infinity of R form a flat 
closed submanifold of M . Thus, M must have a flat closed submanifold of dimension at 
least I. 



Let & C ~# be a closed flat submanifold of maximal dimension. To simplify the visual 
appearance of formulas, denote by F the z-image of J^, and by Fy the orthogonal comple- 
mentary subspace to F in § n at the point r(y). If dim J? = n — 1, then ^ = M and the 
proof is completed. 

Let us turn now to the case of dim < n — 1. Since (^#, i) is locally convex and does not 
have points of strict convexity, there is an open set ^> C ^# that contains & such that 
is of the form <9(0 • F x C F ) fl S n , where C F is a salient (n — 1 — dim ^-dimensional 
convex cone in O-F^. Note that i~ x (F^~ f]i(^)) is a locally convex hypersurface in F- 1 , which 
is strictly convex at y. If dim F -1 = 2, then z is conical - the pointed cone C is the product of 
a linear (n — l)-subspace in IR n+1 and two rays (with origin at 0) in a complimentary linear 
2-subspace. 

Suppose now dimi 7 ^ > 2. Upon applying van Heijenoort's theorem to the map 
c y i : i~ l (Fy fl z'(^#)) — > c^F- 1 ), we see that c y i i~ 1 (F- L fl is a complete convex 

hypersurface in c y (F y L ) C T^. But this is true for all z/ G J^. Furthermore, if y,y' E ^ 
are sufficiently close, the local convexity implies that M y = i i~ l (Fy fl z 



and My/ = 



/(Fy- fl z(*/#)) are isometric: an isometry can be chosen as a minimal rotation p yy i of § n 
in ]R ra+1 around the orthogonal complement to the plane spanned by vectors r(y) and r(y'), 
which maps r(y) to r(y'). Since & is compact, the manifolds M y are isometric for all values of 
y G & . Also, for any z/, z/ G the isometry mapping p TO / from M y to My preserves the dis- 
tance of each point of M y to F. If the map i : z'~ 1 (F y L ni(^)) — > F^ is strictly convex at some 
point z 7^ y, then any geodesic segment I z = i(J? z ), where J? z is a flat open 1-submanifold 
through z G is transversal to Fy. Therefore, any support plane Hj z for through 
I z intersects with a support plane for (^#, i) through F over a proper subset of F. But then 
for any neighborhood ^ y of y there will be y' G ty£ y D such that p^y will move r(z) to a 
point lying on the other side of Hj z with respect to the convex witness of z, resulting in a 
contradiction. Thus, for all y G & the hypersurface c y i iy 1 (F ± fl C c y (F y L ) cannot 

be strictly convex at any z ^ y, which means c y i z'" 1 ^ 1 - fl z(^#)) 
in c y (Fy). This implies that % is conical. ■ 



is a convex salient cone 



The notion of convex part, introduced by Van Heijenoort (1952), happens to be very useful 
in working with local convexity. A convex part of , i) at a point of strict convexity i(o) 
(i.e. containing a point of strict convexity i(o)), is a connected subset C of with z'~ 1 C 

open in ^jf, such that 

(1) i^di^C) — H fl i(^), where F is some hyperplane in X™, not passing through i(o), 

(2) C lies on the boundary of a convex body K c such that dK c C C U F. 
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The set i~ C is called an abstract convex part at o and is denoted by C. We will denote 
i{di~ x C) by reldC and we will call H D Kc the lid of the convex part C. Let H a be a 
supporting hyperplane at i{6). Let us call the open half-space defined by H that contains 
C\i(o) the positive half-space and denote it by . 

Theorem 5 (after van Heijenoort) If an immersion of an (n — 1) -manifold (n > 3) into X™ 
has a point of strict convexity, it has a convex part containing this point. 

Proof, van Heijenoort's proof works for X n , n > 3, without changes. ■ 

Let us fix a supporting hyperplane H Q and consider a family of hyperplanes such that: (1) 
they are pairwise disjoint in Hq and their intersections with Hq form a partition of Hq, (2) 
they are all orthogonal to a line (circle) L, transversal to H , (3) they do not contain i(o). 
Let us call such a family the fiber bundle {H}^ L Ho ) of the positive half-space defined by L 
and H . Denote by A > the length of the segment on the line (circle) L from I = L D H 
to H G {H}(l,h ), called the height of H. The elements of {H}(l,h ) can then be indexed 
by their heights. We use H\ for the hyperplane of height A. For X n = R" the range of A is 
(0, +oo) and for X n = S n it is (0, 7r). If there exists a convex part at i(o), then {H x }(l,h ) 
defines a family {C\} of convex parts at i(o), ordered by inclusion. Each such convex part 
is "squeezed" between H and H\ and inherits the height from its H\. Let us now consider 
the union C Q of all convex parts at i{6): we want to prove that C Q is a convex part itself. 
This statement was proven by van Heijenoort for X n = W 1 (n > 2) and we will prove it for 
S n (n > 3). 

Theorem 6 Let i : ^ — ► § ra be as in Theorem^ and let C = i(G) be a convex part at i(o), 
defined by a hyperplane H\ C S n (X ^ 0) from a fiber bundle {Hx}(l,h ), where L _L H Q . 
Suppose rel<9C is the boundary of a convex set S in H\. Then either S is the i-image of a 
topological disk § fdimS < n — 1) in ^# and = C U §, or C is a proper subset of a larger 
convex part at i(o) and defined by the same bundle {Hx}(l,h )- 

Proof. We will now prove this theorem by a perturbation argument, which reduces the 
spherical case to that of IR n . Since relc^S = rel<9C, S C H\ (A > 0), and dimS" = n — 1, 
we conclude that convC D H Q = i(o). Since, by Lemma [lj M is arcwise connected, all of 
conv C, except for i(o), lies in the positive halfspace . Thus, there is a hyperplane H in § n , 
arbitrarily close to H D and orthogonal to L, such that C lies in an open halfspace H + defined 
by H. Let p be the pole of § n with respect to H that lies in H+, and let c p : H + — ► T p be the 
central projection map on the tangent plane T p C M n+1 . The set i~ 1 (H + n is obviously 
a manifold. The map c p i is a locally-convex immersion of i~ 1 (H + H into T p with a point 
of strict convexity, i(o). Any Cauchy sequence on under the c p i-metric is also a Cauchy 
sequence under the i-metric. Thus the immersion map c p i : i~ 1 (H + H — > T p is complete 
and, therefore, satisfies the conditions of Theorem It maps any (spherical) convex part 
centered at i(o) onto a Euclidean convex part; it also maps the fiber bundle {Hx}(l,h ) to 
a fiber bundle in T p . In the case of X™ = M n the statement of the theorem is known (van 
Heijenoort). We conclude that either = C U S, or C is a proper subset of a larger convex 
part centered at i(o) and defined by the same bundle {H\\^l,h )- ■ 
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4 Prom Local to Global Convexity on the Sphere 

Let us recall that Van Heijenoort proved that a complete locally convex immersion / of a 
connected manifold ^# (dim ^# = n — 1) into M. n (n > 3) is the boundary of a convex body, 
if / has a point of strict convexity. For n = 3 this result, according to van Heijenoort, follows 
from four theorems in Alexandrov's book (1948). Jonker & Norman (1973) proved that if / 
does not have a point of strict convexity, f{^) is the direct affine product of a plane locally 
convex curve and a subspace L = M. n ~ 2 of M. n . 

Theorem 7 Let % : j$ — ► §™ in > 2) &e an immersion of a connected (n — 1) -manifold ^ , 
satisfying the following conditions: 1) ^# is complete with respect to the i-metric, 2) (./#, i) 
is locally convex at all points of \M ' . Then i : ^# — ► S n is an embedding onto the boundary 
of a convex body in § n . 

Proof. If points of strict convexity are absent, then by Theorem [5] the map i : j& — > S n 
is an embedding onto the boundary of a convex body. Thus, we can assume that a point 
of strict convexity exists. If o G ^ is such a point, then by Theorem [5] there is a convex 
part containing i(o). Consider the union C a of all convex parts at i(o). Denote by C the least 
upper bound of the heights of convex parts defined by {H x }(l,h ) at i(o). By Theorem [6] 
dC a C H^ni(^) and rel9C is the boundary of a closed convex set D in H^. Two mutually 
excluding cases are possible. 

Case 1 : dimD < n — 1. Then, following the argument of van Heijenoort (Part 2: pp. 239- 
230, Part 5: p. 241, and Part 3: II on p. 231), we conclude that C Q U D is the homeomorphic 
i-image of the (n — l)-sphere i~ l (C ) U i~ 1 (D) C where i~ 1 (C Q ) is a maximal abstract 
convex part at o. Since is connected, i~ 1 (C Q ) Ui^ 1 (D) = and i : ^# — > § n is a convex 
embedding of 

Case 2 : dim D = n — 1. By Theorem[6] C Q is either a proper subset of a larger convex part, or 
C , together with the lid D, is the homeomorphic i-image of ^ . Since the former alternative 
is excluded by definition of C Q , i : ^ — > § n is a convex embedding onto <9(CU D). ■ 

Note that the statement of the above theworem is invalid for n = 2. For example, although 
the 1-surface in § 2 depicted in Figure 1 is locally convex at all points, it does not bound any 
convex set on S> 2 . 



5 Locally convex PL-surfaces 

Let 7 be a fixed star-finite semi-regular cell-partition of j$ . Recall that in our terminology a 
cell is always homeomorphic to an open ball. We say that r is locally convex at a cell C G T 
if it is locally convex at each point of C. 

Theorem 8 Let r : ^ — > IR n {n > 2) be a complete cellwise-flat PL-realization of a con- 
nected manifold ~# (dim^ = n — 1) such that 
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Fig. 1. Locally-convex PL-hypersurface in § 2 , which is not convex 

1) r is locally convex in at least one point of each (n — 3) -cell. 

2) r is bounded or strictly locally convex in at least one point of \M ' . 

Then r : ^# — > M. n is an embedding onto the boundary of a convex body defined by (possibly 
infinitely many) affine inequalities. 

Proof. Upon invoking the definition of cellwise-flat PL-realization and that of local convexity, 
we conclude that r is an immersion. We know that r : ^# — > M. n is locally convex at all 
(n — 3)-cells. Since r is a cellwise-flat PL-realization, r is also locally- convex at all rf-cells 
for d > n — 3. If r is bound ed, then convr(^#) is bounded and compact. By Straszewicz's 
theorem ( Rockafellar . 1990), the set convr(^#) has an exposed point. As r is complete, this 



exposed point must the r image of some x G j$ . Since r is locally convex at x, it is also 
strictly locally convex at x. If we prove that r is locally convex at all cells, by Theorem [3] 
the map r : j$ — ► MJ 1 is a convex embedding. We proceed by reverse induction in cell's 
dimension. Suppose we have shown that r : ^# — > K n is locally convex at each fc-cell, where 
< k < n — 3 (and therefore at all cells of higher dimensions). If n — 3 = 0, the proof is 
finished. So, let n > 4 and let us consider a (k — l)-cell F G 7. Consider r(StarF) D §p, 
where is a sufficiently small (n — fc)-sphere lying in an affine subspace complementary to 
aff r(F) and centered at some point of r(F). Also note that dimSp = n — k > 2. The map 
r : ^# — > M. n is locally convex at F if and only if the hypersurface D r(StarF) in is 
convex. Since r : ^ — > MJ 1 is locally convex at each fc-cell, the surface D StarF is locally 
convex at each vertex and therefore locally convex everywhere. The set r _1 (Spflr(Star F)) is 
complete in the r-metric and thus, by Theorem [7J the surface §p R r(Star F) is an embedded 
convex hypersurface in §p. Notice that the condition dim§^ = n — k > 2 is essential to the 
applicability of Theorem [7] (see Figure 1 for a locally convex surface in § 2 which is not a 
convex surface in S 2 ). So, r : ^# — * M. n is locally convex at F. 

To recap, the above induction argument shows that r is locally convex at all vertices, and, 
therefore, at all points. The metric induced by r is indeed complete. Upon applying Van 
Heijenoort's Theorem [3] and Theorem [71 we conclude that r is an embedding and that r(^f) 
is the boundary of a convex body. ■ 
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Corollary 9 Let r : ^ — > IR n in > 2) be a complete cellwise-flat PL-realization of a 
connected manifold ^ (dim^ = n — 1). Suppose r is bounded or is strictly locally convex 
in at least one point. If \M ', IP) is (n — 3) -simple, i.e. exactly three (n — l)-ce//s mafe contact 
at each {n — 3) -cell, then r{./&) is the boundary of a convex polyhedron. 



6 Convexity Checker for PL-hypersurfaces 

In this section we present a polynomial-time algorithm for checking the convexity of any PL- 
realization r : ^# — > W 1 (n > 3) of a semi-regular cell-partition T of a connected compact 
(n — l)-manifold The map r under testing is assumed to be cellwise-flat (see Section [2]), 
which implies that each cell C G IP is homeomorphicly mapped by r to an open subset of an 
afline subspace of dimension dim C. We do not assume that the realization is an immersion: 
if it is not an immersion, the algorithm will detect this. We do not make any generic position 
assumptions. 

In describing the algorithm we assume that certain combinatorial and geometric information 
is readily available. This input information is exactly what should be kept by a convex 
hull computer if it is to use our verification procedure. Later we discuss the complexity 
of extracting the necessary input information from PL-surface descriptions given in some 
typical formats. 

If any of the subprocedures return 1 or "false", the main procedure returns "false" as the 
final answer. The idea of the algorithm is to check that the immersion and the local convexity 
properties hold at each corner. For each corner C this check is reduced, roughly speaking, 
to the verification of convexity of a certain cone K(r, StarC) in r(C) ± = IR 3 , which is 
constructed from the poset of Star C and the restriction of r : — > M™ to Star C. Such a 
cone is not unique - for example, any non-singular afline transformation of K(r, Star C) is 
just as good as K(r, Star C). This reduction from the star of a corner to a cone in R 3 is done 
by the procedure Reduce-to-3D. 



6.1 Input Conventions 

Let X be a subset of W 1 . We use aff X to denote the afline subspace spanned by X and 
> I 

aff X to denote the linear subspace {x — x' x,x' G aff X}. For a (not necessarily convex) 

> 

polytope IIcK™ with a face F an inner normal at F is any vector n in aff II such that for 
any point p on F there is e > such that p + en G int II. The normal n is called a Euclidean 
normal if n _L aff F. 

As a reminder, we refer to the (n — 3)- and (n — 2)-cells of the partition as corners and ridges 
respectively, and we refer to (n — l)-faces as facets. We also use these terms to refer to the 
realizations of these cells in M n . Mathematically, the input is given as follows: 
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(1) the subposet CP[n — 3, n — 2, n — 1] C CP of corners, ridges, and facets where it is known 
in advance which are which; 

(2) a Euclidean inner normal to R at C for each ridge-corner incidence (R, C); 

(3) a Euclidean inner normal to F at C for each facet-corner incidence (F, C). 

The data in (1) will be referred to as combinatorial, and that, described in (2) and (3), 
as linear-algebraic. We assume that each vector in the linear-algebraic data "knows" the 
corresponding abstract cells in CP , and that each abstract cell in CP[n — 3,n — 2,n — 1] 
"knows" all normal vectors related to it. The input data-structure can be implemented as 
a double-linked adjacency list, with appropriate attribute fields for dimensional and linear- 
algebraic data. Namely, we can create an adjacency list for the directed (multi-) graph whose 
vertex set consists of elements of CP[n — 3, n — 2, n — 1] and whose edge set consists of all 
ordered pairs (C,C) and (C',C), where C -< C or C -< C in 3>[n - 3, n - 2, n - 1]. When the 
input is available in this form we say that the input is given in the standard form. 

In applications a PL-hypersurface is usually specified by a subposet of the face poset, which 
includes the vertices or the facets or both; it is normally equipped either with the coordinates 
of vertices or with the equations (or unequalities) for the facets. Suppose now the input is 
given as the poset CP[0,n — 3, n — 2, n — 1], equipped with the coordinates of the vertices; 
in this case we will say that the input is given in traditional form. If the partition CP is a 
triangulation, then the linear-algebraic data required for our algorithm (standard form) can 
be produced in linear time in f n _3„_2, which is also polynomial in the total bit size of the 
input. More generally, if the face numbers of facets of , CP) are bounded by a universal 
constant (in f n _3„_ 2 ), the linear-algebraic data for the standard form of the input can be 
computed by using 0(f n _ 2n _ 3 ) field arithmetic operations. 

6. 2 Preprocessing 

By preprocessing in the context of problems of verification of geometric properties we mean 
any computation that does not depend on the geometric realization (in our case r) , but only 
on the topology or combinatorics of the object (in our case - the pair (^#, CP)). 

Since jtft is a manifold, the facets of (^#, CP) making contact at a corner are "glued" to each 
other in a circular fashion. Same can be said about the ridges. The circular structure of the 
stars of (n — 3)-cells implies that for each (n — 3)-cell C we have f n „2(StarC) = f n _3(StarC). 
The last formula implies that for the whole IP we have f n _3 n _2 = f n _3„_i. More properly, 
a topologist would say that the "links" (defined via the 1-skeleton of the dual partition, 
a well-known construction going back to H. Poincare: see Seifert & Threlfall, 1980) of the 
corners are circles. These circles can be thought of as polygons whose vertices correspond to 
the facets of 7 and edges to the ridges of CP. To apply our algorithm we need to determine 
a cyclic order of ridges around each corner. Such an order is unique up to the choice of 
direction.. To apply our algorithm for different realizations of the same cell-partition of ^# it 
is reasonable to maintain a cyclic order of ridges around each corner; this can be considered 
as preprocessing. 
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6.3 Algorithm 

The main procedure Convexity-Checker is given in Algorithm 1. Convexity-Checker works on 
a stack Corners, in which we put all (n — 3)-faces of (^#, T) prior to starting. Subroutines 
Reduce-to-3D and Check-if-Cone-Convex are used by Convexity-Checker, and subroutine Is- 
Folded is used by Reduce-to-3D. Before describing the working of the algorithm, let us define 
a couple of auxiliary notions. 

For any m > 2 a graph with vertex set V = {0, 1, . . . , m} (m > 2) and edge set E = 
{(01), . . . , (0m)} U {(12), . . . , (m — 1 m), (m, 1)} is called the m-wheel graph and is denoted 
by W m . Vertices RimV(W m ) = {l,...,n} are called rim vertices; edges RimE(W m ) = 
{(12), . . . ,(m—lm), (m, 1)} are called rim ridges. Edges {(01), . . . , (0m)} are called spokes. 
Vertex is called the center of W m . Let p : V — > R 3 be a realization of the vertex set of W rn 
in R 3 . If m > 3, then we can assign to each 3-cycle (Oii + 1) (where i + 1 is taken mod m) in 
W m a geometric simplex in R 3 with the vertices p(0), p{i), and p(i + 1). Therefore, the map 
p : V — > R 3 produces a simplicial surface with boundary. With a slight abuse of terminology 
we will say that a realization p : V — > R 3 of the wheel graph W m (V,E) (where m > 3) is 
convex if p is injective and the resulting simplicial 2-surface, which we denote by p[W m ], is 
convex at p(0). While we may encounter Euclidean realizations of the vertex set of W 2 , we 
will not have a need to associate a surface in R 3 with such realizations. 

Once a corner C is popped from the stack, a pair (W m , n) = (W m [C], n[C]) is created. This 
pair consists of the wheel graph W m = W m [C], which encodes the combinatorics of StarC, 
and an array of vectors n[C], whose elements are the inner normals to the r-realizations of 
ridges and facets of Star C at the face r(C). Namely, m is the number of ridges meeting at C, 
the center of W m corresponds to C, the rim vertices correspond to the ridges of Star C, the 
rim edges correspond to the facets of Star C, and the "spokes" correspond to the corner-ridge 
incidences. 



Algorithm 1 Convexity-Checker 
while Corners ^ do 

Pop an (n — 3)-face C from Corners 
Create (W m ,n) for StarC 

if Reduce-to-3D(VF m , n) = 1 then return false > not convex 

else 

if Reduce-to-3D(W m ,n) ^ then 
(W m ,p) <- Reduce-to-3D(VF m ,n) 

if Check-if-Cone-Convex-in-3D(W / m ,p)=/a/se then return false > not convex 
end if 
end if 
end if 
end while 

return true >yes, convex 
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6.4 Reduce-to-3D and its Subroutines 



For a corner C of T) let {0, 1, ... , m} be the vertices of the wheel graph W m = W m [C], 
with as the center of W m (which corresponds to the conrer C), and let [1, . . . , to] be a cycli- 
cally ordered list of the rim vertices of W m (which correspond to the ridges of Star C) . The in- 
put to Reduce-to-3D consists of: the wheel graph W m ; the Euclidean inner normals ri\, n rn 
for the r-images of ridges of Star C at r(C); the Euclidean inner normals n 12 , n ml for the 
r-images of the facets of StarC at r{C). The arrays of normal vectors can also be thought 
of as a map n : RimE{W m ) U RimV(W m ) -> W n . 

Reduce-to-3D uses a subroutine Is-Folded, which checks for violations of the immersion as- 
sumption in cases where all normals are in the same plane. The output of Reduce-to-3D is 
either one of {0, 1} or a realization of W m in R 3 , where the center of W m is mapped to the 
origin. Output 1 means that the input (W m , n) is inconsistent with our assumptions about 
the integrity of the input data or the immersion assumptions. Output means that Star C 
has passed the local convexity check. 

6.4.I Is-Folded 

Let v, u, w be three coplanar non-zero vectors. The ordered triple (v, u, w) defines a plane 
angle at the origin in the following way: v and w span the two extreme rays of the angle, 
while u is an interior vector of the angle - i.e., the function of u is to specify which of the two 
open subsets defined by v and w is interior to the angle. We denote such angle by (v\u\w). 
Note that (v\u\w) = (w\u\v). 

Is-Folded takes as input a 5-tuple of coplanar non-zero vectors (a,b,c,d,e), where a,b,c 
are pairwise distinct and c,d,e are pairwise distinct. Is-Folded returns true if the interiors 
of angles {a\b\c) and (c|d|e) overlap and false otherwise. For example, Figure 2 shows the 
case where "folding" takes place: angle (c|d|e) "folds over" the angle (a\b\c). 




b 



Fig. 2. Angles (a\b\c) and (c\d\e) overlap. 

In the pseudocode of this procedure we will use a boolean predicate P(v\ui,u 2 \w), which 
is defined for any 4-tuple of coplanar vectors v,ui,u 2 ,w, where v and w are distinct and 
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Ui ^ v, Ui ^ w for i — 1,2. P(v|iti, u 2 |iu) is /aZse if (u|iti|iw) 7^ (t>|u 2 |u>) (Figure 3, right) 
and true otherwise (Figure 3, left). 




Fig. 3. Left: P(v\ui,U2\w) = false. Right: P(v\ui,U2\w) = true. 



Input: a, b, c,d,eE W 1 , where dimspanja, b, c, d, e} = 2, \{a, b,c}\ = 3, |{c, d,e}\ = 3 
Output: boolean 

Algorithm Is-Folded 

if a and e define the same ray then 

if P(c\b, d\a) = true then return true 

else return false 

end if 
end if 

if P(a\b, e\c) = false and P(c\d, e\a) = true and P(b\c, d\e) = true then return false 
else return true 
end if 



The following algorithm shows how to compute P(v\ui,U2\w) via standard linear algebra. 
For any ordered pair of vectors [ei,e 2 ], such that {ei,e 2 } C span{v, u±, w 2 , w}, we use 
sgn[ei,e 2 ] to denote the orientation of [ei,e 2 ] with respect to some fixed orientation of 
span{t>, Ui, u 2 , w}). 
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Algorithm P(v\ui,U2\w) 
s <— sgn[v, w] 

if sgn[iti,to] = sgn[v,iti] = s then 

if sgn[u 2 ,w] = sga.[v,u 2 ] = s then return true 

else return false 

end if 
else 

if sgn[n, 2 , w] = sgn[v, u 2 ] = s then return false 
else return true 
end if 
end if 

6.4-2 Procedure Reduce-to-3D 

Input: W m : m-wheel graph, n : RimE(W m ) URimV(W m ) — > M. n . In the pseudocode we use 
rii for n(i) and n ii+ i for n(i i + 

Output: one of {0; 1; (W m ,p) where p : V(W m ) -> R 3 } ) 
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Algorithm 2 Reduce-to-3D 



e\ <— ni, S <— {ei} // S is a maximal independent set of vectors in (aff r(C)) 1 " 
if rank{ni,n 2 } = 1 then 

if n 2 = Xrii for some A > then return 1 

end if 

else e 2 <— ra 2 , S <— S U {e 2 } 
if m = 2 then 

if {rankjnx, n 12 , n 21 } = 2} and {sgnfr*!, n 12 ] = sgn[ni, n 2 i]} then return 1 
else return 
end if 
end if 
end if 

% <— 3 // rim- vertex counter 

while i ^ 1 mod m and |S| < 3 do 

if rank{ei, e 2 , n^} < 2 then / /n*, e l5 e 2 are all in one plane 

if nj = Ae! or = Ae 2 for some A > then return 1 

end if 

if |S| = 1 then / /e\ and e 2 are collinear and contraoriented 
if sgn[e 2 ,ni 2 ] = sgn[e 2 ,n 2m ] then return 1 
else 

if m = 3 then return 
else e 2 <— n 3 , S <— S U e 2 
end if 
end if 

else //in this case we know |S| = 2 

if ls-Folded(e 1 , n 12 , e 2 , nj_! j, rij) = true then return 1 
end if 
end if 
else e 3 <— n u S <— {e 3 } 
end if 
i <- i + 1 
end while 

if i = 1 mod m or |S| < 3 then return 
end if 

p(0)<- (0,0,0) 
for j = 1 to m do 

<- (ei • nj, e 2 • nj, e 3 • n^) 

end for 
return (W m ,p) 
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6.5 Procedure Check-if-Cone-Convex 



Denote by p[W m ] the simplicial surface that results from the map p : V(W m ) — > M 3 (see 
Subsection 6.3). Informally speaking, we test the surface p[W m ] for convexity by going around 
the wheel and checking whether the following two conditions are satisfied or not. The first 
condition will be the conical PL-surface pfW 7 ™] "turning in the same direction" each time 
we increment the index: here, the formal meaning of "the same direction" is captured by 
the linear-algebraic notion of orientation of a frame. Recall that the sign of a list [vi,v 2 , 1*3] 
of three vectors in R 3 is the sign of the determinant of the 3x3 matrix whose i-th row 
is Vi, which we denote by sgn[t?i, V2, V3]. If Vi,V2,Vs G M™ with n > 3, then to define the 
sign of the triple [v±, V2,Vs] we need to fix an orientation in a 3-subspace of IR n containing 
{vi,V2,Vs}. Let sgn[p(l),p(2),p(3)] = s ^ 0. Formally, the first condition is that for every 
i G V(W m ) the sign sgn\p(i),p(i + l),p(i + 2)] must be s or 0; the zero sign corresponds to the 
case where the surface "continues straight", i.e., vectors p(i),p(i + l),p(i + 2) lie in one plane. 
The second condition will be p[W m } not intersecting itself. This can be captured by checking 
the sign of every triple of the form [p(l) , p(i) , p(i + 1)] is or s = sgn[p(l),p(2),p(3)]. Thus, 
roughly speaking, the first condition ensures that the convexity is not lost due to a turn in 
the wrong direction, while the second condition guarantees the cone will not intersect itself. 
Finally, the correctness of Check-if-Cone-Convex hinges on the following lemma, whose proof 
we omit. 

Lemma 10 Let p : V{W m ) — > R 3 be a realization of the m-wheel graph (m > 3) which does 
not map any two consecutive spokes into the same ray. Suppose sgn[p(l),p(2),p(3)] 7^ 0. 
Then p[W m ] lies on the boundary of a convex cone if and only if for all < i < m — 1: 



The input to Check-if-Cone-Convex-in-3D is wheel graph W m (m > 1), each of whose vertices 
v is equipped with a corresponding point p{y) in R 3 . For notational simplicity we assume 
that the vertex set V(W m ) of W m is {0, 1 ... , m}, where is the center of W m , and [1, . . . m] 
is a cyclic order on the rim vertices. In addition, we assume that vertex is realized at the 
origin 0. 



Input: W m : wheel graph, p : V(W m ) — > IR 3 



Output: boolean 



1) sgn[p(i),p(i + l),p(i + 2)] 



sgn[p(l),p(2),p(3)] and 



2)sgn\p(l),p(i + l),p(i + 2)\ 



sgn[p(l),p(2),p(3)]. 
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Algorithm 3 Check-if-Cone-Convex 

sign <- sgn[ei,e 2 ,e 3 ] 
for j from 2 to m do 

if {sgn[p(j),p(j + l),p(j + 2)] = -sign} or {sgn[ei, e 2 ,p(j + 2 )] = -sitfra} then 
return false 

end if 
end for 
return true 



7 Complexity Analysis 

In this section f& denotes the number of fc-cells of (^#, D 5 ), and the number of incidences 
between fc-cells and /-cells of (^#, T); if 5" is a subposet of D 5 , then ffc(J") denotes the number 
of /c-cells in J". Note that all linear algebra in the algorithm is essentially reduced to com- 
parisons of signs of lists of at most three n-vectors; we will refer to any such calculation as 
a sign computation. Unless mentioned otherwise, as e.g. in the next Subsection, we assume 
that the input is in the standard form. 

(1) Building the wheel graph for StarC takes time linear in the number of ridges of Star C. 

(2) Since Corners is accessed at most f n _ 3 times, Reduce-to-3D is called at most f n _ 3 times. 

(3) Reduce-to-3D(IU m , n) requires at most 0(m) sign computations. 

(4) Check-if-Cone-Convex requires at most 0(m) sign computations. 

(5) Is- Folded requires a constant number of sign computations. 

1 Suppose the algorithm uses the field arithmetic (+, — , x,+) and each arithmetic operation 
has unit cost. This model is realistic when real computations are conducted with floating 
point arithmetic. If n is fixed, the complexity of the algorithm is 0(f„_ 3n _2) = 0(f n _ 3n _i). 
To estimate the complexity in the case where n is one of the parameters describing the 
input size, we need to estimate the contributions of sign computations in (3)-(5). Notice 
that any sign computation in (3)- (5) deals with, at most, six n-vectors. Since standard 
linear-algebraic procedures over a field can be used, the complexity of the algorithm is 
0(nf n _ 3n _ 2 ). 

2 What follows is a discussion of the complexity in the cases where no floating point error 
can be tolerated. Let R be the base ring of the computational model: i.e., all numerical 
input data (such as the coordinates of vertices, the coefficients of normals to (n — l)-faces 
etc.) come from R. Furthermore, we assume that Z < R < M. in the partial order of rings. 
When we discuss the degrees of the polynomial predicates evaluated by the algorithm, 
we consider them as polynomials with integer coefficients in the input parameters. In this 
context the phrase arithmetic operation stands for any ring-theoretic operation (+, — , x). 

Case 1: the dimension n is fixed. In this case all linear- algebraic computations can 
be done via determinants. Using determinants has an advantage of keeping the degrees 
of evaluated polynomial predicates at bay. Moreover, since in our algorithm the largest 
determinants are 3x3, the highest degree of evaluated predicates is 3. Thus, the arithmetic 



21 



complexity of the algorithm is 0(f ra _3„_2) and the algorithm evaluates at most 0(f n _3„_2) 
polynomial predicates of degree 3. 

Case 2: the dimension n is not fixed. If n is not too large, the linear-algebraic compu- 
tations can still be done via direct determinant evaluations. In each computation we are 
dealing with at most three n- vectors, which means that we may have to evaluate ( \) 



3x3 determinants to find a minor of maximal rank. Thus, the total arithmetic complexity 
of the algorithm is 0(n 3 f„_ 3n _ 2 ). The case of large values of n is considered in the next 
paragraph. 

1.1 Exact Computations over Z. 

In here we consider the case of exact computation. The dimension n is not fixed and R — Z. 
We are now interested in the bit complexity of our algorithm, e.g., in the multitape Turing 
machine model. Note that the same techniques can be used for sign computations in our 
algorithm. Since each sign computation involves no more than 6 vectors, the bit complexity 
of each sign computation using Yap's (2002) ramification of the Bachem-Kannan algorithm 
is (9(nMfe(L)), where Mb(x) is the bit-complexity of multiplication of two integers of binary 
sizes not exceeding x and L is a bound on the binary size of the coefficients of the vectors (see 
(Yap, 2002) for details). Then the total complexity of the algorithm is 0(nf n _3„_ 2 M{,(L)). 
Devillers et al. (1998) have shown that any convexity checker, whose work does not depend 
on the nature of R, has to evaluate at least one polynomial of degree n - however, this lower 
bound is mandatory only for those checkers that work the same way for any R. 

Theorem 11 Let r : jj( — > W l be a cellwise-flat PL-realization of a manifold (^C, IP) of di- 
mension n—1. Suppose the input is in the standard form and all normals have integer coor- 
dinates of binary size not exceeding L. There exists a polynomial time algorithm for checking 
convexity ofr : — > M, n with (multitape Turing machine) complexity of 0(nf n _ 3 n _ 2 M fe (L)) . 

Now, let us consider the situation where the input is given in the traditional form, i.e. as the 
poset CP[0, n — 3, n — 2, n — 1], equipped with the coordinates of the vertices. If we have no 
restrictions on the combinatorics and geometry of the geometric realization of cells of (IP, ^#) , 
then it is very difficult, or even impossible, to construct Euclidean inner normals from given 
data. Let us assume that the partition 7 is simplicial. In order to do sign computations, we 
need first write down Euclidean inner normals for all corner-facet and corner-ridge incidences. 
For each such incidence we have to deal with roughly n vectors of length n. Computing a 
Euclidean normal is then reduced to a fidning a non-zero solution for a homogeneous system 
Mx = where M is at most n by n matrix. We can use Yap's (2002) version of the 
Bachem-Kannan algorithm to compute (upper triangular) Hermit's Normal Form for the 
system Mx = 0. Then a non-zero solution vector of at most polynomial size can be found 
in polynomial time by using standard techniques of linear algebra: we just work our way 
from the bottom of the normalized matrix up until all %i S £1X6 found. Alternatively, one can 
use a polynomial algorithm in Yap (2002: Sec. 10.8-10.9), based on repeated application of 
Bachem-Kannan algorithm, to further reduce the system to Smith Normal Form and then 
find a solution. Furthermore, to reduce the complexity, we can deal with each corner C in 
the following way. If v , — ,v n are the vertices of r(C), then we first find Hermit's normal 
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form for the matrix [v\ — Vq, ...,v n — Vq] and then, for each P, where P is a ridge or a 
facet incident to C, compute an integral Euclidean normal vector to r(P) at r(C). Then the 
complexity of all linear-algebraic computations for StarC is dominated by the complexity 
of finding Hermit's normal form for the matrix [vi — v , ...,v n — v ], which is 0(n 3 Mb(L)) 
(Yap, 2002). Thus , the total complexity is 0(n 3 f„_ 3 „_ 2 M 6 (L)). 

If n = 3, then the corners are the vertices and the required normals are easy to produce. 
Now, suppose n > 3. How large the coefficients of the integral Euclidean normals, discussed 
above, can be? It is obviously possible to produce each such Euclidean normal as a vector 
whose coordinates are polyno mials of deg ree at most n — 3 in the coordinates of the vertices. 



Furthermore, Siegel (see e.g. lYapl . 120021 ) proved that a homogeneous system of k linear 
equations with n variables over Z has a non-zero solution where each component is bounded 

k 

in absolute value by 1 + (nA)^-^ (for us k = n — 3) where A is the largest of the absolute 
values of the coefficients. Siegel also showed this bound could not be improved. When n is 
small enough, a vector satisfying Siegel's bound can be found by classical methods of lattice 
reduction (no efficient methods for finding such a vector are known for large n). If normals 
satsifying Siegel's bound are used in the algorithmic precedures given above, then the largest 
integers that may appear in sign computations via determinants are of the order \(nA ) n ~ 3 
- where A ^ 1 is a constant slightly larger than 1, and A is twice the largest of the absolute 
values of the vertex coordinates. 



Theorem 12 Let r : ^# — > M n be a cellwise-flat PL-realization of a simplicial manifold 
(^#, y) of dimension n — 1. Suppose the input consists of the poset CP[0, n — 3, n — 2, n — 1] 
equipped with the coordinates of the vertices, and that for each vertex v = r(v) of r(^#) we 
have v e Z™ and \v\ < 2 L . There exists a polynomial time algorithm for checking convexity 
of r : ^ — > IR n with (multitape Turing machine) complexity of 0(n 3 f n _3„_2Mb(L)). 

The input requirements in the above theorem can be relaxed. If we know only "J" = !P[0, n — 
3, n — 2] (or J" = CP[0, n — 3, n — 1]) together with the circular order of facets (or ridges) at 
all ridges, then 7[0,n — 3,n — 2,n — 1] can be computed at no extra cost. 



7.2 Surfaces in M 3 



The algorithm runs in linear time in the number of vertices when ^ is spherical. However, 
a sequence of non-spherical PL-manifolds of dimension 2 can have the edge number growing 
quadratically in f . Thus, it is desirable to check the topological type of the input by just 
counting 1-cells (edges) in CP[0, 1]: once their number exceeds cfo, where c is some constant 
which is easy to calculate, we stop and declare the input non-convex. This check helps 
preserve the O(fo) running time bound for PL-surfaces in IR 3 . One may wonder if such a 
check is necessary, as it seems very likely our algorithm will quickly encounter a non-convex 
vertex, if the input surface is homeomorphic to a sphere-with-handles or sphere-with-Mobuis- 
strips. Surprisingly, Betke & Gritzmann (1984), proved that any orientable non-spherical 
connected closed 2-manifold can be PL-embedded into M 3 so that it has exactly 5 non- 
convex vertices but no fewer! The problem of determining the minimal possible number of 
non-convex vertices in a PL-immersion of a non-orientable closed 2-manifold is open. 
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In the case of M 3 the requirements on the combinatorial part of the input can be somewhat 
relaxed: in what follows we show it is sufficient to know only CP[0, 1], which is the 1-skeleton 
graph of CP). First, the planarity of this graph can be checked in O(f ) time (Lempel 
et al, 1967). For a planar graph we can also determine the faces in linear time - i.e., in 
O(fo) time we can create the face-nodes, where each face-node is double-linked to its edge- 
nodes. Once we know the faces in terms of their edges, we can double link each face-node 
to the vertices-nodes of all of its edges. Because of the sphericity of (^#, CP) the latter task 
takes O(f ) time. Thus, the adjacency list representing CP[0, 1, 2] can be constructed from the 
adjacency list representing CP[0, 1] in O(f ) time. 

The case of M 3 is a rather special one. First, n = 3 is the smallest dimension for which the 
techniques of this paper apply. Second, even in the case of M. n the convexity test for each 
corner is reduced to testing convexity of a section of the star of this corner, which is essentially 
equivalent to testing convexity of a cone in IR 3 . Another important consideration is that in 
application a 2-surface in IR 3 is normally specified by its combinatorics and the coordinates of 
the vertices or equations for the facets: it is therefore important to specify how our algorithm 
can be applied when the imput is given in the traditional form. Namely, suppose we are given 
CP[0, 1,2] equipped with the coordinates of the vertices v 1 = r(v i), . . . , v fo = r(v fo ). The 
corner- ridge normals are then just vectors Vi — Vj. The question remains how to find corner- 
facet normals, i.e. vectors pointing from the vertices of the facets into the interiors of the 
facets. This is easy if it is known that the facets are convex. Otherwise we have the following 
algorithmic problem. Let Ck be the /c-cycle graph. Consider a rectilinear embedding r of 
C n in an affine plane A C M 3 - the pair (Ck,r) defines a 2-dimensional polytope P(Ck,r) 
whose boundary is r(Ck) (here Ck is regarded as PL-manifold). Let v be a vertex of Ck- 
The problem is to find a non-zero vector n e A such that r{y) + en lies in the interior 
of P(Ck,r). This problem can be solved in time O(k); solving this problem for all facets 
will require 0(fo) ring-arithemtic operations. Thus, there is no difference in time-complexity 
betwen the standard and traditional forms of the input for n — 3. We will now restate the 
observations made in this Section in the following theorem. 

Theorem 13 Let r : ^# — > IR 3 be a cellwise-flat PL-realization of a 2-manifold (^#, CP). 
Suppose we are given the 1-skeleton of , CP) equipped with the coordinates of the vertices. 
Let L be the upper bound on the bit sizes of the coordinates of the vertices. There exists an 
algorithm for checking convexity of r : ^ — > IR n with (multitape Turing machine) complexity 
of O(f M 6 (L)). 

8 Conclusions 

This paper describes a local approach to convexity verification of PL-hypersurfaces. The main 
theoretical result of the paper is a characterization of global convexity of a PL-hypersurface 
in IR n in terms of the local convexity properties of the surface at its (n — 3)-faces. Building on 
this approach we give a polynomial-time convexity checking algorithm that can be applied 
for any closed hypersurface. The approach presented in this paper can be generalized to 
piecewise-polynomial surfaces of small degree. 
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